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ABSTRACT 



We show that there are no dynamic screening corrections to the Salpeter's 
enhancement factor in the weak-screening hmit. 

Subject headings: nuclear reactions 
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1. Introduction 

Continuing observations of solar neutrinos and helioseismological analyses of the solar 
interior (Bahcall 1989 and references therein) are becoming more accurate and warrant 
accurate theoretical expressions for thermonuclear reaction rates in astrophysical plasmas. 
As shown by Salpeter (1954), decreased electrostatic repulsion between reacting ions caused 
by the Debye-Hiickel screening leads to a noticeable increase in reaction rates. The reaction 
rate enhancement factor is 

w = exp( ^^ ), (1) 

where Ry) is the Debye radius. The above expression is valid in the weak-screening limit 
which is defined by TRd ^ Z^Z^e^ ■ 

Carraro, Schafer & Koonin (1988) noted that if the Gamow energy of the reaction is 
high enough, reacting nuclei may have velocities much higher than thermal ion velocities. 
The plasma ions response will be suppressed, because at high enough velocities the ion 
part of the plasma dielectric constant is nearly zero ("dynamic screening"). The authors 
suggested that in such cases, only electron part of the screening should be included into (1), 
and Eq. (1) should read 

At even higher Gamow energies, the electron part of the dielectric constant also goes to 
zero, and, according to this argument, the screening enhancement effect disappears. 

We will show that although dynamic screening can be important for certain laboratory 
experiments, the nearly precise thermodynamic equilibrium of stellar plasmas guarantees 
that Salpeter's expression, Eq. (1), is valid independently of the Gamow energy (§3). Before 
that, in §2, we summarize the formal part of the argument of Carraro et al. (1988). 
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2. What is Dynamic Screening? 

Consider a test charge Zic moving through a plasma with velocity v. The plasma 
response to the test particle is described by the dielectric constant e which is a function of 
both the wavenumber, /c, and the frequency, u. For the purpose of reaction rate calculation, 
we need to know the electrostatic potential, 00, created by the plasma at the test particle 
location. The enhancement factor is then (Salpeter, 1954) 

w = exp{-Z2e(f)o/T), (3) 

for a Zi — Z2 reaction. Now, knowing e one calculates 0o as 

For t> = 0, we have e = 1 + {kRo)"'^, and (4) gives 0o = —Zie/Rc,. When inserted into Eq. 
(3), this result gives the Salpeter enhancement factor, Eq. (1). 

If reacting ions are fast moving, with velocities v in the interval VTi < v < vtb, the 
dielectric constant in Eq. (4) is close to the static dielectric constant of the electron part 
only, that is e ^ 1 + (A;i?£))~^/2. With this dielectric constant, Eq. (4) gives a \/2 smaller 
potential (f)Q, which translates into a smaller enhancement factor, Eq. (2). 

The above analysis seems well based, but is in contradiction with the following simple 
argument. The Gibbs probability distribution, p, for the plasma is 

p~exp(-/3E^-/3E— )> (5) 

where (3 is the inverse temperature. This distribution is factorable, 

pr^eM-pj:^)x^M-pj:—)- (6) 

That is, distributions in velocity and configuration spaces are decoupled: fast moving 
charges are just as screened as others. Salpeter's result, Eq. (1), can be derived starting 
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from the probability distribution for spatial configurations, and thus does not depend on 
the Gamow energy. 

3. A Hypothetical Reaction 

The strongest disagreement between Salpeter's expression and the dynamic screening 
prediction would occur in an extreme case of very high Gamow energy. Then interacting 
nuclei are moving at very high velocities, and the plasma does not respond at all. We will 
show that even in this case Salpeter's expression is correct. 

The plasma does not respond to fast moving nuclei Zx and Z2, but the nuclei move in a 
thermal (fluctuating) electrostatic field, 0(r). We will show that reaction rates in a random 
potential are increased, and for thermal fluctuations this effect is correctly accounted for by 
the Salpeter expression. 

In thermal equilibrium, the local densities of the nuclei are 

«i,2(r) = Ci,2exp(-/3Zi,2e0(r)). (7) 

The local reaction rate is i^ni(r)n2(r), where i^ is a multiplicative factor that depends only 
on temperature. The average reaction rate is, for a gaussian random field 0, 

R = K<nx{r)n2{r)> = i^CiC2exp(^/3V(Zi + ^2)' < 0^ >), (8) 

where we have taken < >= with no loss of generality. Now, the average densities of the 
reacting nuclei are related to the normalization constants C\^2 by 

<ni,2> = Ci,2exp(i/3VZ22<02>), (9) 

and, in terms of mean densities, Eq. (8) reads 

R = K <ni><n2> exp{f3'^e'^ZiZ2 < 0^ >). (10) 
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This corresponds to an enhancement factor 

w = exp{p'^e'^ZiZ2 < 0^ >). (11) 

Equation (11) is valid for an arbitrary gaussian random potential. To make a 
connection with Salpeter's formula, we have to calculate the size of the thermal electrostatic 
fluctuations, < 0^ >. The fluctuating electric field, E, is given by the equipartition law, 

where the notation < E'^ >k is defined by 

The potential fluctuation is < 0^ >fc=< E'^ >k /k'^, and Eq. (12) gives 

< 02 >^^ (47rT)A;-2(l + pRl)-\ (14) 

Now, 



and Eq. (11) gives 



w = exp( ^^ ). (16) 



This is exactly Salpeter's result. It was obtained without explicitly introducing the concept 
of screening. 
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